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We present a simple reaction kinetics model to describe the polymer synthesis used by Lusignan 
et al. 0] to produce randomly branched polymers in the vulcanization class. Numerical solution of 
the rate equations gives probabilities for different connections in the final product, which we use to 
generate a numerical ensemble of representative molecules. All structural quantities probed in ref. Q 
are in quantitative agreement with our results for the entire range of molecular weights considered. 
However, with detailed topological information available in our calculations, our estimate of the 
'rheologically relevant' linear segment length is smaller than that estimated in ref [lj. We use a 
numerical method based on tube model of polymer melts to calculate the rheological properties 
of such molecules. Results are in good agreement with experiment, except that in the case of the 
largest molecular weight samples our estimate of the zero-shear viscosity is significantly lower than 
the experimental findings. Using acid concentration as an indicator for closeness to the gelation 
transition, we show that the high-molecular-weight polymers considered are at the limit of mean- 
field behavior - which possibly is the reason for this disagreement. For a truly mean-field gelation 
class of model polymers, we numerically calculate the rheological properties for a range of segment 
lengths. Our calculations show that the tube theory with dynamical dilation predicts that, very 
close to the gelation limit, contribution to viscosity for this class of polymers is dominated by 
the contribution from constraint-release Rouse motion and the final viscosity exponent approaches 
Rouse- like value. 



PACS numbers: 82.70.Gg, 83.10.Kn, 02.70.-c 



I. INTRODUCTION 

Polycondensation reactions that generate branched 
polymers lead to progressively larger molecules as a func- 
tion of the conversion. The same thing happens during 
chemical cross linking (vulcanization). At a critical ex- 
tent of the reaction or density of bonds p c , the size of 
the largest molecule spans the system and this is termed 
as the gel point 0,^0. Close to the gel point, static 
properties of the system exhibit a scaling form. Defining 
e _ Ip~p<=I > the number fraction $(Af) of the molar mass 
M falls off as a power law, 



$(M) ~ M- T f(M/M char ), 



(1) 



where / is a cut-off function and the characteristic mass 
M char diverges as 



Mcha 
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Different moments of the molar mass distribution also 
diverge as the gel point is approached. In particular, the 
weight-averaged molar mass Mw diverges as 



Mi 



w 



(3) 



The static exponents r, a and 7 depend on the univer- 
sality class for a given system. When the molecules in 
the melt overlap strongly (as in vulcanization of long lin- 
ear molecules), the exponents belong to the mean- field 
universality class and are described by Flory-Stockmayer 
theory 0,0]. In this case, the exponents can be calcu- 
lated analytically with r = 5/2, a = 1/2 and 7 = 1. 



Polymerization of small multifunctional groups lead to 
the critical percolation gelation class. Though exponents 
for this class cannot be calculated analytically, good es- 
timates for the exponents are known from simulations 


Close to the gel point, under some circumstances, the 
rhe olog ical properties also obey scaling forms [9. Hoi HH 
IT2I ITaTIT^. Il5j . The shear relaxation modulus G(t) is a 
power law in time (t) and the complex viscosity T)*(u>) is 
a power law in frequency (w): 
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t~ 



and 



(4) 



The zero-shear viscosity 77 diverges with exponent s and 
the recoverable compliance J® diverges with exponent t: 
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and 



(5) 



The dynamic exponents are not derivable from the 
static ones without further assumptions |l6j and for en- 
tangled polymers the effective exponents (since the relax- 
ation is only approximately a power law \L7\ ) depend on 
the length of the linear segments between branch-points. 
By considering simple rules for relaxation as a function 
of the seniority variables of the segments in a branched 
material, the variation of the dynamic exponents as a 
function of segment length has been estimated 01 ■ A de- 
tailed calculation which takes care of the molecular topol- 
ogy without such approximations is missing. In a recent 
publication Q , we used a numerical method to calculate 
the relaxation of arbitrarily branched material within the 



broad framework of tube theor y |l 8l an d its extensions 
to handle constraint release [lj, |2iJ, an d constraint- 
release Rouse motion [zj H3 • Relaxation of branch-on- 
branch architectures were included in a manner which 
respects the polydispersity both in length and in topol- 
ogy. In this paper we attempt to use such a numerical 
scheme to calculate the rheological relaxation function 
and the dynamic exponents close to the gel point. 

The rest of the paper is organized as follows. In 
a recent paper, Lusignan el al. Q reported synthesis 
and characterization of a series of randomly branched 
polyesters which are in the mean-field gelation class. Af- 
ter a brief description of their reaction scheme fsec. ITT|l . 
we use a simple kinetic model to determine the vari- 
ous probabilities for the connectivity in the final prod- 
uct. Using such probabilities we generate representative 
molecules and characterize the static structural proper- 
ties and compare them with the experimental findings 
(sec. Illlf) . We provide a brief qualitative description of 
the numerical method of 2] to calculate the rheological 
properties of branched entangled polymers in sec. IIV Al 
In sec. IIV Hi we calculate the linear rheological response 
of the molecular ensembles and compare with the exper- 
imental results. The average inter-branch-point segment 
length in these polymers depends on the extent of esterifi- 
cation and the estimate of this length is subject to errors. 
In sec. we consider a simplified ensemble of molecules 
which have predetermined average segment length and 
investigate the segment length dependence of the dynam- 
ical exponents. We conclude the paper by recapitulating 
the main findings of this study and stressing the questions 
this work raises on our understanding of the relaxation 
of highly branched polymers at the longest timescales. 



II. KINETIC MODELING FOR BRANCHED 
PTMG POLYMERS 
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FIG. 1: (a)-(c) Molecular structure of the reactants in poly- 
condensation reaction considered in Qj (d) Schematic descrip- 
tion of the esterification involved in the synthesis. 



In a recent study, Lusignan el al. [l| considered a 
polycondcnsation reaction of a polytetramethylene gly- 
col (PTMG) oligomer with number-averaged molar mass 
M N = 2900 g/mol, trimethylolpropane (TMP) and 
adipic acid (AD). The two OH groups at the ends 
of PTMG and three OH groups at the ends of TMP 
molecules (see figure ^) react with the two acid groups 
of AD. Thus AD works as a bridging molecule connect- 
ing the TMP and PTMG molecules. The trifunctional- 
ity of TMP molecules leads to branching. FASCAT 4100 
(monobutyl tin oxide) was used as catalyst which be- 
comes incorporated in the final product. In our simplified 
description, we assume that the catalyst simply increases 
the reaction rate without changing the final product. 

PTMG and TMP molecules were mixed in a 3:1 molar 
ratio. The fraction of AD was controlled to generate sam- 
ples of different molecular weights. The acid numbers at 
the end of the reaction were near zero, - signifying com- 
plete conversion. The intrinsic viscosity of the molecules 
so generated show a transition from the linear-like behav- 
ior at the low molecular weights to randomly branched 
behavior at the high molecular weights. The crossover 
between this region was found at Mx = 66000 g/mol. 
The high value of Mx compared to mass of the oligomers 
(PTMG) indicates that a large fraction of the TMP 
molecules has one of the OH groups unreacted. 

In our modeling, we consider that all AD molecules 
present react completely with some OH group. As in the 
synthesis, we consider molar ratios n PTMG : ?i T mp = 3:1. 
To keep the number of free parameters to a minimum, we 
assume the same rate constants (ki) of esterification for 
the OH groups on PTMG oligomers and that on the unre- 
acted TMP molecules. Once one of the OH groups on the 
TMP molecule has reacted, it can affect the rate constant 
for the second OH group because of the small separation 
of the OH groups on the TMP molecules. Thus we as- 
sume different rate constants (£3) for rate constants 
of esterification of OH group on TMP provided one (two) 
of the OH groups has already reacted. In what follows, 
we denote unreacted TMP molecules as TMPo and TMP 
molecules with n of their OH groups reacted as TMP„. 
The kinetic equations considered are 
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; 1 and [ 

^cqoh] refer to the concentrations of 
unreacted OH groups on the PTMG molecules and num- 
ber of unreacted COOH groups on the acid respectively. 



2 



[n" MP ] refers to the concentration of TMP m species. We 
fix the acid concentration as n AD = / a n PTMG , where f a 
determines the extent of stoichiometric imbalance and 
hence the extent of the reaction. The rate equations 
are solved numerically with the initial condition t = 0, 
Pcooh] = 2n AD , [n^ MP ] = n TMP , [^o™ G ] = 2tt.p T mg 9-nd 
[ n ™p] = f° r m being 1,2 or 3. At the end of the 
reaction [n COOH ] = 0. Thus scaling the concentrations 
by [n COOH ], the procedure involves numerical integration 
from [n COOH ] = 1 to [n COOH ] = 0. 



and the probability weight of each individual molecule is 
simply the inverse of the total number of molecules so 
generated. For the rheological response, the molecules 
contribute with this probability weight. We have here 
assumed that there are no ring molecules and that the 
reactivities are independent of the size of the molecule - 
both of which assumptions are expected to break down 
as the gel point is approached. Also our analysis depends 
on the assumption of spatial homogeneity (continuously 
stirred reaction scheme). 
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FIG. 2: Fraction of TMP molecules having 1, 2 or 3 of the OH 
groups reacted as a function of average molar mass, kz/ki is 
assumed to be 0.8. 

To reduce the number of free parameters further, we 
assume that the presence of one reacted OH group on a 
TMP molecule lowers the reaction rate by a certain frac- 
tion and the presence of two reacted OH groups inhibit 
the reaction rate of the third OH group independently 
(hi, = /c|). For a given acid concentration and value 
of k2/k\, the numerical solution of the rate equations 
yields the probabilities of having different reacted species 
in the final product (fig. |2J). From these probabilities, we 
generate an ensemble of representative molecules by first 
selecting a species (PTMG, TMP„ or AD) with proba- 
bilities given by their respective weight fractions. Any 
unreacted acid group reacts with an OH group on either 
a PTMG molecule or a TMP„ molecule with the proba- 
bilities from the solution of eqn.[|J] For TMP„ molecules, 
n of the end groups are attached to AD molecules. For 
PTMG molecules, ends are attached to an acid group 
to have the probability of reacted OH groups on PTMG 
the same as that given by the solution of eqn. |SJ Since 
the initial species is selected on a weight basis, in this 
way we generate a weight-biased molecular distribution 



III. STATIC STRUCTURE OF BRANCHED 
PTMG POLYMERS 
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FIG. 3: Determining M c h ar '■ the line represents the best 
fit (f){M) ~ M~ :i/2 exp(-M/M char ) indicating the exponent 
t = 5/2 and M char = 2 x 10 5 
in this plot. 



g/mol. k^jkx is fixed to be 0.8 



For a given choice of fe/fci, the acid concentration f a is 
varied to generate a series of different A% ensembles. For 
mean-field gelation ensemble, the cut-off function in the 
molar mass distribution f(M/M c har) in eq.^is explicitly 
known to be |5( 



f(M/M char ) = exp 



M 



(J) 



Using this cut-off function and assuming that the ex- 
ponent r = 5/2 in eqn. ^ we use a two-parameter fit 
to determine M c ^ ar from the tail region of the molar 
mass distribution (fig. Note that because our en- 
semble is generated on weight basis, we fit a function 



$(m) - Ar 3/2 



exp 



M 



2M C , 



corresponding to r = 5/2. 



The small-mass end of the distribution (not shown in the 



3 



figure) does not conform to this form and shows noisy 
features due to the finite size of the oligomers used dur- 
ing synthesis. 
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FIG. 4: M c h a r as a function of Mw for several choices of 
fe/fci ratios superposed with the experimental data. 

In fig. 0] we show M c h ar as a function of Mw for three 
different choices of k-z/kx- For &2 = k\, the simulation 
values of M c h ar are consistently higher than the experi- 
mental points, while for It2 = 0.7fci, the simulation values 
are consistently lower. For hi = 0.8fci, the simulation re- 
sults closely match the experimental values. 

Ref. measured the intrinsic flow viscosity, [rj\, and 
found a crossover from linear-like behavior ([77] ~ M - 8 ) 
at low molar mass to randomly branched behavior ( [77] ~ 
M 0A5 ) at high molar mass. When all the samples of 
different Mw were considered together, the crossover 
of these two behaviors was found at Mx — 6.6 x 10 4 . 
Without knowledge of the detailed interaction among the 
monomers, it is not possible to compute the intrinsic vis- 
cosity. The intrinsic viscosity should depend linearly with 
the radius of gyration in a good solvent - which again is 
beyond our ability to calculate. However, it is reasonable 
to assume that the radius of gyration in a good solvent 
will be directly related to that in the 8 solution, and 
in particular that the crossover molar mass between the 
linear-like and the branched scaling will be the same. 
From the numerical ensemble of the polymers, we used 
Kramers theorem to calculate this ideal radius of gy- 
ration. 

From the numerical ensemble of the molecules, we form 
histograms of molecules with respect to the mass of the 
molecules. For each bin in the histograms we calculate 
the average radius of gyration. For these calculations, 
we assume that the Kuhn mass is 74 g/mol and the 
results are in units of Kuhn length b. In fig. [S] we plot the 
radius of gyration for two different molar mass samples 
(symbols), both generated with fe/fei = 0.8. Also shown 
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FIG. 5: Ideal radius of gyration for &2 = 0.8fci and Mw = 
15200 (filled circles) and M w = 226000 g/mol (open trian- 
gles). The line shows the combined behavior of all the data 
sets from samples with different Mw together, b is the Kuhn 
length and Kuhn mass is assumed to be 74 g/mol. 



is the radius of gyration when all the different molar mass 
samples are considered together (line). The individual 
Mw ensembles roughly fall on this merged distribution 
line. A difference shows up only at the high molar mass 
limit - where the lower Mw ensembles do not have any 
entries. 

In fig. El we plot the radius of gyration when all sam- 
ples are considered at fc 2 = 0.8fci. At the low-molar- 
mass end, the data fits the form R 2 g ~ M corresponding 
to linear polymers. The high-molar-mass end fits the 
form R 2 ~ M 1 / 2 corresponding to randomly branched 
polymers. The crossover of these two behaviors is found 
by extrapolating the fits at Mx — 65800 g/mol. In- 
creasing (decreasing) the value of ki leads to lowering 
(raising) Mx- For comparison, h% = 1 and = 0.7 
respectively corresponds to Mx — 47500 g/mol and 
Mx = 89000 g/mol. Since the same value of fe/fci fits 
both this crossover behavior and the variation of M c ^ ar 
with Mw with the experimental results, only results with 
this value of fo/fci are shown in the rest of the paper. 
The quantitative agreement with quite different exper- 
imental results using the same parameter suggests that 
our simplistic assumptions about the reaction kinetics are 
probably close to reality. 

From the cross-over in radius of gyration, one might 
conclude that the typical linear segment length is about 
66000 g/mol (indeed, this is the conclusion in p| from the 
intrinsic viscosity data). With the detailed topological 
connectivity at our disposal, we can probe the segment 
length between branch points in a different way. A linear 
segment can be made of PTMG oligomers connected by 
acid groups alone - or with intervening TMP molecules 
with only two of the three OH groups reacted. Such 
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FIG. 6: Crossover from linear to randomly branched behavior 
from radius of gyration. The symbols correspond to the radius 
of gyration determined by making a histogram in mass from 
all the different Mw samples considered together. The solid 
line and the broken line corresponds to slopes 1 (linear) and 
1/2 (randomly branched) respectively. The arrow indicates 
crossover at 65800 g/mol. fe/fci is fixed at 0.8. 
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FIG. 7: Probability distribution of linear segments at a fixed 
acid concentration. p(Mn,s) decays exponentially at large 
Mn,s- The line denotes an exponential fit in the shown range. 



TMP connectors will have small side-arms (compared to 
entangled molecular mass) which will still behave as lin- 
ear segments in rheological measurements. We therefore 
add the mass of such small side-arms to the backbone. 
This does not remove segments which are too small to be 
rheologically important but still are connected at all ends 
(for example, an H molecule with the cross bar formed 



by two TMP molecules will behave just like a four arm 
star). For this reason we take an alternative route than 
a simple average of the masses of the arms (we believe 
that this alternative route is also Theologically most rel- 
evant). Any random association to form linear segments 
attains a Flory distribution (most probable distribution) 
at mass scales much larger than the constituent elements. 
The probability of having a segment of length M a is given 
by p(M a ) = c exp(— M/Mjv,s)j with c being an constant 
and, Mn,s being the number-averaged molar mass of the 
segments. From a histogram of linear segments, we deter- 
mine p(M a and using an exponential fit at large M a de- 
termine the number-averaged molar mass Mjv,s (fig. 0). 
The weight-averaged molar mass for Flory distribution is 
twice Mn,s- 
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FIG. 8: Estimate of number averaged molar mass of linear 
segments between branch points as a function of Mw- The 
error bars are estimated from standard error in the exponen- 
tial fit of p{M a ) (fig. El. 

In fig [S] we plot Mm,s as a function of Mw- The error 
bars correspond to the error estimates in exponential fit 
of M(a) (fig. 0). determining the addition probability. 
At large Mw, the estimate of the linear segment length 
from this approach is ~ 7200 g/mol, which is almost an 
order smaller than Mx determined from the crossover of 
radius of gyration. This is due to the fact that lightly 
branched material like stars or combs, which dominate 
the mid-range in the mass distribution, have a radius of 
gyration which is closer to that of linear polymers with 
the same molar mass than to that of randomly branched 
polymers. 

As an estimate for closeness to the gelation transition, 
we define 
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(8) 



where, is the critical acid concentration where the 
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FIG. 9: Dependence of characteristic mass M c h aT on the acid 
fraction /„. The line is a linear fit for high Mw samples. The 
intersection point with the zero x-axis gives the critical acid 
concentration f£, where the characteristic mass diverges. 



characteristic mass diverges. Close to the gel point, the 



characteristic mass M c j lar scales as M c h ar ~ e . Thus 
the plot of 1/yMchar versus f a shows a linear behavior 
for large f a (fig. EJ • The point at which the line crosses 
the zero x-axis (infinite M c i iar ) determines = 1.173(1). 
For k-i/ki = 0.8, the sample with average molar mass 
220 Kg/mol corresponds to f a = 1.156 giving e ~ 0.014. 
The size of the largest branched molecule provides a char- 
acteristic length £ and mean-field theory provides a self- 
consistency test by requiring that the molecules of this 
characteristic size should overlap sufficiently. In 3 dimen- 
sions, this leads to a critical value of the extent of the 
reaction e c ~ iV -1 / 3 below which the largest molecules 
no longer overlap significantly and the exponents change 
from the mean-field results Q, 0]. Taking the linear 
segment length -Mjv,s = 7200 g/mol and Kuhn mass to 
be 74 g/mol, there are on average approximately 100 
Kuhn segments between branch points. This estimate 
gives e c ~ 100 -1 / 3 ~ 0.2, so that the highest molar mass 
samples are much closer to the gelation transition than 
the critical value and are therefore expected to show non- 
mean-field behavior. We will meet the rheological conse- 
quences of this critical behavior in the following. 



IV. RHEOLOGICAL RESPONSE OF 
BRANCHED PTMG POLYMERS 

A. Computational rheology 

To estimate the rheological behavior of the entangled 
branched PTMG polymers, we employ a numerical ap- 
proach dm. 

For details, the reader is referred to [2j. 



We summarize the procedure qualitatively here for com- 
pleteness. The numerical approach is based on tube the- 
ory 0, which replaces the topological entanglements 
from neighboring chains by a hypothetical tube surround- 
ing a given chain. After a small strain, the stress is re- 
laxed by the escape of the chains from the old tube con- 
straints. This connects the stress relaxation to the sur- 
vival probability of the chains in their respective tubes. 
In a polymer melt, since all the polymers are in motion, 
the tube constraint itself is not fixed over time. This 
constraint release is handled by the dynamic dilation hy- 
pothesis 0, l2fl l2l| , which postulates a simple relation 
between the tube diameter and the amount of unrelaxed 
material. A free end monomer relaxes part of its tube 
constraint at short times by constraint-release Rouse mo- 
tion. At later times, the entropic potential, which itself 
evolves due to constraint release, leads to a first-passage 
time approach |26j| . The contribution from a collapsed 
side arm is modeled by including increased friction on 
the backbone, estimated from the time of collapse and 
the current tube diameter as a length-scale for diffusive 
hops from an Einstein relation. For branch-on-branch ar- 
chitectures, the relaxation leads to a multi-dimensional 
Kramers' first-passage problem. We simplify this by re- 
casting it to an effective one-dimensional problem which 
has the required Rouse scaling, respects topological con- 
nectivity and gives correct result at some special known 
limits [3 • A linear or effectively linear (branched material 
with collapsed side arms) chain can relax by reptation. 
When a large amount of material relaxes quickly, such 
that the dynamically dilated tube increases in diameter 
faster than the rate permitted by Rouse relaxation, the 
effective orientational constraint responds more slowly by 
constraint release Rouse motion and the dynamic dila- 
tion is modified [22I l23l |. In addition, we include contri- 
butions from the Rouse motion inside the tube and fast 
forced redistribution of material at the early stages of the 
relaxation |27j . 

In computational rheology, starting from a numerical 
ensemble of molecules, tube survival probability in dis- 
crete (logarithmic) time is followed after an imaginary 
step strain [2j|25l|28|,|2!|. At each of these time steps, the 
amount of unrelaxed material 4>t and the effective amount 
of tube constraint 4>st is stored. Since the visco-elastic 
polymers have a very broad spectrum of relaxation, we 
assume that the amount of material relaxed in each time 
step contribute as independent modes in the stress relax- 
ation modulus G(t). Thus after all of the molecules have 
relaxed completely, G(t) is calculated as a sum over all 
these independent modes. The complex modulus G*(oj) 
at frequency lj is defined by 



G» 



ijj I G(t)exp(—iu)t)dt. 
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(9) 



The real and imaginary parts of G*{uj), storage modulus 
G (a;) and the dissipative modulus G (w) respectively, 
are of particular interest since they are measured in os- 
cillatory shear experiments. The zero-shear viscosity r\ is 
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calculated from 



B. Dynamic exponents for branched PTMG 



770 = irm 



CO 

and the steady-state compliance J° is calculated from 



(10) 



J° = lim 
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(11) 



All the integrations are replaced by sums over discrete 
timesteps of the relaxation. 

The calculations have a few free parameters. The 
material-dependent parameter of entanglement molar 
mass M e is related to the plateau modulus Go by j^lj 



_ ApRT 
Go= 5 MT' 



(12) 



where, p is the polymer density and T is the temperature. 
The timescale is set by the entanglement time r e which 
is the Rouse time of the chain segment between entangle- 
ments. When the molar mass of the segments is scaled by 
M e and the time is scaled by r e , in the approximation of 
tube theory, polymers of the same topology but of differ- 
ent chemical composition relax the same way. We assume 
that, for a side-arm relaxing completely at certain time 
t a , at times much larger than t a , the motion of the asso- 
ciated branch-point can be modeled as a simple diffusion 
process with hop size pa at the timescale of t a . Here a 
is the tube diameter and p is a numerical factor. We use 
p 2 = 1/40 as used in to fit a wide range of different 
experimental data. The dynamic dilation hypothesis as- 
sumes that the effective tube diameter depends on the 
amount of unrelaxed material and the effective number 
of entanglements associated with a segment of length Z 
scales as Z — > <f)^ a . We choose the dynamic dilation 
exponent a = 1 in our calculations. 

For the class of polymers considered in this study, the 
number of branches on a given molecule can be quite 
large. Also a large number of molecules need to be con- 
sidered to ensure that a single massive molecule does 
not affect the results disproportionately. In our ap- 
proximations, the relaxation of the different molecules 
are coupled only via the amount of unrelaxed mate- 
rial (fit- This enables us to divide the ensemble of 
molecules in several subsystems and follow the relax- 
ation process independently at each step, communicat- 
ing the local <pt to other processors at the end of each 
time step. The minimal communication needed makes 
the parallel-code scale almost perfectly with the num- 
ber of processors and most importantly allows us to 
probe closer to the gelation transition, where the mem- 
ory requirement becomes larger than available on single 
processors. The source code, precompiled executables 
and documentation of the program are available from 
http : //sourcef orge .net/projects/bob-rheology. 



To calculate the dynamic properties of the branched 
PTMG molecules, we generated ensembles of 5 x 10 s 
molecules at each Mw considered and followed the re- 
laxation after a small step-strain. In the absence of high- 
frequency measurements for this material, we take M e 
and T e as free parameters - fitted to describe the dynamic 
properties. Without the complications of occasional ester 
groups from the esterification and the butyl side groups 
from the TMP molecules, the present polymers resem- 
ble polytetrahydrofuran (PTHF). For PTHF, treating it 
as an alternating copolymer of ethylene and ethylene ox- 
ide, the estimate of the entangled molecular weight is 
M e ~ 1420 g/mol We use this value as our rough 
first guess for M e and fix r e by matching the zero-shear 
viscosity with experimental results at the intermediate 
molar mass range of the experimental data. 




FIG. 10: Zero-shear viscosity i/asa function of weight- average 
molar mass. Filled circles are the experimental data from 
p|. The triangles and the squares are the results from our 
calculations using M e = 1420 and 1080 g/mol respectively. 
The error bars for the 1080 g/mol calculations represent the 
effect of uncertainty in r\ corresponding to the experimental 
uncertainty of 10% in determining Mw- The line is a power- 
law fit to the M e — 1080 g/mol results in the intermediate 
mass region with the exponent 3.95. 

In fig. ^|we plot the zero shear viscosity rj for differ- 
ent values of Mw- The filled circles are the experimental 
results from lj. The triangles are results from our cal- 
culations with M e — 1420 g/mol and r e = 1.8e — 7 s. 
For this choice of M e , the viscosity increases slowly with 
Mw compared to the experimental data. The squares 
are results with M e = 1080 g/mol and r e = 3.45e — 8 s. 
This choice of M e is able to reproduce the experimental 
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viscosity data over half a decade in My/. A power-law 
fit in the intermediate range (shown as a dashed line in 
the figure) gives the viscosity exponent s — 3.95(2). In 
the experiments, there is an uncertainty of 10% in de- 
termining Mw- The error bars in the M e = 1080 g/mol 
data show the associated uncertainty in viscosity. At the 
largest Myy, our calculations and the experimental data 
show opposite trends. The viscosity from our calcula- 
tions shows a trend of lowering of the exponent at the 
largest Mw, while the experimental data shows a sharp 
increase. For rest of the results in this section, we use 
M e = 1080 g/mol. 
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FIG. 11: Complex viscosity t)*(uj) for selected Mw- Lines 
are from calculations in this study, symbols are experimental 
data from [l|. 

Fig. ^2 shows the frequency dependence of the com- 
plex viscosity r]*(uj) for several different values of Mw- 
Symbols represent experimental data from |lj in the in- 
termediate mass range, where zero-shear viscosity from 
our calculation matches with the experimental values. 
rj*{uj) shows an approximate power-law behavior with 
exponent 1 — u. Away from the gelation transition, this 
power-law window is limited. We fitted power laws in 
the frequency range 10-100 s _1 to estimate the exponent. 
Since 1/Mw ~ e, we plotted u{Mw) from such fitting as 
a function of 1/Mw (fig- I12fl . Linear extrapolation to 
1/Mw — > gives the limiting value of u = 0.305(1), 
which corresponds closely to the experimental value of 
u = 0.31(2). 



V. MEAN-FIELD GELATION ENSEMBLE 

The segment length for branched PTMG considered in 
the earlier part of the paper is largely determined by the 
size of the oligomer used to synthesize the polymers. In 
this section we turn to a hypothetical series of polymers 
which fall in the category of mean-field gelation class. 
We consider linear molecules of type A which are Flory 
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FIG. 12: Variation of the relaxation exponent u with Mw- 
The line shows the extrapolation used for large Mw to find 
the value of u at the gel transition. 



distributed and tri-functional groups B which have zero 
mass. The molecules are formed by the rule that A reacts 
with B and neither A- A nor B-B reactions are allowed. 
Furthermore, we assume that, at the end of the reaction, 
all B bonds are attached to some A molecules. Thus, as 
in the case of branched PTMG polymers, extent of the re- 
action is determined by stoichiometric mismatch. As be- 
fore, we assume that there are no closed loops. The final 
distribution of the molecules are described by only two 
parameters: Mjv.s, number-averaged molar mass of the 
linear segments and pb, the branching probability. The 
molecules are generated by selecting the first strand with 
a Flory distribution of length Mm,s and adding Flory 
distributed branches recursively on both ends with prob- 
ability pi,. 

The static properties of these molecules can be solved 
analytically. The characteristic molar mass diverges 
when pt is 0.5 (e = (0.5 — pb)/0.5). Using seniority 
variables to approximately describe the hierarchical re- 
laxation, ref. |32l ] found that the entangled contribution 
to the terminal relaxation time of this class of polymers 
does not diverge at the percolation threshold. Their 
calculation did not include the constraint release Rouse 
modes, contributions from which will still be divergent 
in the absence of a diverging entangled contribution. 
For the calculations presented in this section, we assume 
M e = 1120 g/mol and r e = 1.05 x 10 -8 s, corresponding 
to high-density polyethylene at 150° C 0- For each value 
of pb and Mjv,s considered, we generate an ensemble of 
2 x 10 5 molecules and follow the relaxation after a step 
strain. To estimate statistical errors involved in our cal- 
culations, for each case we repeat the calculation 3 times 
with different sets of molecules (generated by different 
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random seeds). 




FIG. 13: Radius of gyration for gelation ensemble with 
Mn,s ~ 22400 g/mol with pb = 0.49. The cross-over mass 
Mx = 314400 is roughly 14 times larger than Mn,s- 

In fig. we plot the mass dependence of the radius 
of gyration for pi, = 0.490 (e = 0.02) and segment length 
Mn,s = 22400 g/mol (number of entanglements between 
branch points N/N e = 20). As in the case of branched 
PTMG, the extrapolated cross-over mass Mx in radius of 
gyration from linear to the randomly branched behavior 
is much larger than Mm,s- Because the segments are 
Flory distributed with out a lower cutoff, the difference 
is even larger in this case (Mx /Mm,s ~ 14). 

Fig. ^^illustrates the procedure followed for estimating 
the apparent relaxation exponent u (when the relaxation 
dynamics are entangled there is no reason to expect true 
power-law behavior, but an apparent power law can hold 
as a good approximation for a sizeable range of relaxation 
timescales 17]). The left subpanel shows the variation of 
the complex viscosity rf with frequency for three dif- 
ferent e for N/N e = 20. At the lowest frequencies, for 
e = 0.4, the terminal relaxation leads to significant de- 
viation from the power-law behavior. For smaller values 
of e, this deviation shifts to smaller frequencies. For dif- 
ferent e, we fit a power law with exponent 1 — u in the 
frequency range 10 2 — 10 4 s _1 . In the right subpanel of 
fig. 1141 we plot the dependence of such apparent u with 
e. For e < 0.1, the values of u shows a linear dependence 
on e. A linear fit was used to estimate the extrapolated 
value of u at e = 0. 

Fig. H31 shows the extrapolated values of u (circles) at 
e = from the procedure outlined in fig. [21 as a func- 
tion of number of entanglements between branch points 
N/N e . The error estimates for u are smaller than the 
size of the symbols. An approximate calculation of hier- 
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FIG. 14: Estimating the exponent u : the left panel shows 
t]*(lj) for N/Ne = 20 and e = 0.4, 0.2, and 0.02. Data for 
different e are fitted separately to a power law with exponent 
1 — i/ in the frequency range 10 2 — 10 4 s _1 . The exponent u 
so obtained, as a function of e, is plotted in the right panel. 
A linear fit for e < 0.1 was used to find the limiting value of 
u at this N/N e . 
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FIG. 15: Variation of (apparent) relaxation exponent u with 
N/N e . The circles are the data obtained from our calcula- 
tions. The dotted line is the prediction of hierarchical relax- 
ation model of and the dashed line is the phenomenologi- 
cal form used in Q . The triangle corresponds to the exponent 
u for the branched PTMG polymers, when the segment length 
is estimated by fitting a Flory distribution to the probability 
distribution of segment lengths. 
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archical relaxation in entangled mean-field gelation tube 
model at the gel point 01 predicted a form of u as 



ib— for N > N e 
N ~ 



(13) 



with ip being a constant. Both theory and experiments 
for N < N ei where unentangled Rouse dynamics domi- 
nate, suggest ip ~ 0.67. The dotted line shows the pre- 
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diction from eqn. ^] Ref. Q uses an empirical function, 



0.67 



3 + 21n(iV/7V e ) 



N < 2N e 



N > 2N e 



(14) 



to describe the dependence of u on N/N e from experi- 
mental data. The dashed line in fig. ^| shows this phe- 
nomenological function. Results from our calculations 
fall roughly midway between the prediction of the ap- 
proximate model [17J and the phenomenological fit to 
data in pi. The significant deviation from eq.^Jis mostly 
because uses Mx as an indicator for Mx,s (so over- 
estimating it) and to some extent because u changes ap- 
preciably as the limit e — *■ is considered ffig. H4f> . When 
plotted against our estimate of linear segment length 
(shown as triangle in fig. H4|) . the exponent u correspond- 
ing to the experiments reported in ref. [l| matches closely 
with our calculations on gelation ensemble. 
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FIG. 17: Viscosity exponent s as a function of N/N e . The 
circles are results from our calculations and the dashed line is 
the phenomenological form of Lusignan et al. 
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FIG. 16: Zero-shear viscosity 77 and recoverable compliance 
Je as a function of e for different N/N e . The left subplot 
shows 77 for N/Ne = 3, 5, 7, 10, and 20. The right subplot 
shows J° for only iV/7V e > 7. 

The left subpanel of fig. ED shows the zero-shear vis- 
cosity 77 as a function of e for different value of N/N e 
(larger N/N e data have higher viscosity at the same e). 
The right subpanel shows the recoverable compliance J° 
for N/N e = 7, 10, 15, and 20. The data shows a large 
amount of scatter. J® can be expressed as the first mo- 
ment of G{t) : J e ° = (J™tG(t)dt) /rf. Thus, J e ° is par- 
ticularly susceptible to the long-time decay of G(t). The 
longest relaxation time is dominated by just a few of the 
high molar mass molecules in our ensemble. Thus the 
variation of J° with the particular ensemble considered 
is large. For small N/N e , the relative contribution from 
this tail region of molar mass distribution is even higher. 
For N/N e < 7, the scatter becomes larger than the value 
of J° and they are neither shown in the figure nor con- 
sidered for further analysis. 



At the smallest values of e plotted in the log-log plot 
in fig^l in double log plot, the slope of 77 with e starts to 
decrease. To probe at even smaller e would require much 
larger computations than used in this study. Instead we 
focus our attention in the range of e between 0.006 and 
0.2, where the viscosity for all values of N/N e shows ap- 
proximate power-law dependence on e. Fig. 1 171 shows the 
viscosity exponent s as a function of N/N e . Also shown 
is the phenomenological form of as dashed line 



1.33 N < 2N e 

2\n(N/N e ) N>2N e . 



(15) 



Results from our calculations show a much sharper in- 
crease of s with N/N e than predicted by this functional 
form. 

In fig. El we plot the recoverable compliance exponent 
t (circles) as a function of N/N e . Because of large scatter 
in J° , the error estimates in this case are large (error bars 
are estimates of error from the variance obtained from 
three independent sets of calculations). Since both 77 and 
J° can be expressed as integrals over G(t), the exponents 
u, s and t are not independent. If G{t) behaves like 
t~ u till the longest relaxation time, one gets a dynamic 
hyperscaling relationship among the exponents 



t 



s + t 



(16) 



The estimates of t from estimates of u and s using this 
hyperscaling relation is shown as the squares in fig. 1181 In 
the range of N/N e , where we have direct estimates of t, 
estimates from the hyperscaling relationship falls below 
the direct estimate. 

To explore why this is so, in fig. El we plot the decay 
of G(t) for different values of e and N/N e — 3. Also 
shown is the limiting power law decay suggested from 
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VI. CONCLUSIONS 
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FIG. 18: Recoverable compliance exponent t as a function of 
N/N e . The circles represent results from a direct power-law 
fit of the form J e ° ~ <T l . The error bars are estimated from 
the variance in J° from three independent sets of molecular 
ensemble. The squares are calculated by using the hyperscal- 
ing relationship (cq. H6t . For N/N e — 3 and 5, the scatter in 
is too large for a direct estimate of t. 
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FIG. 19: G(t) (symbols) for N/N e = 3 for e = 0.1, 0.08, 
0.06, 0.04 and 0.02. As e is made smaller, G(t) decays slowly. 
The vertical dotted lines show the range in which the com- 
plex viscosity is fitted (in frequency) to find the value of the 
exponent u. The solid line represents this limiting power law 
decay G(t) ~ t~ u with the indicated slope u = 0.4253. 



fitting the complex viscosity data. For even the lowest 
e studied here, the power-law behavior holds in only a 
small window of the relaxation time and the contribution 
in i] or J° from decay at times larger than this power-law 
window is not negligible. Thus, the dynamic hyperscaling 
relationship holds only approximately (fig. I18fl . 



We have presented a simple kinetic modeling scheme 
for the gelation ensemble polymer synthesis in [l|. With 
just one global fitting parameter describing the branch- 
ing chemistry, we are able quantitatively to reproduce 
the variation of characteristic molar maSS Mchar &s & 
function of Mw and the behavior of intrinsic viscosity 
as a function of molar mass. With the detailed knowl- 
edge of the molecular topology in our calculations, our 
estimate of the average segment length between branch- 
points is much lower than estimated in |l|. We have 
used a numerical technique based on the tube theory of 
polymer melts to calculate the dynamic response of the 
polymers in the linear response regime. For intermedi- 
ate ranges of Mw, both the complex viscosity if(u>) and 
the zero-shear viscosity rj matches with the experimental 
findings. For the largest Mw considered, r\ in our calcu- 
lations is significantly lower than the experimental data. 
At those Mw, our estimate of closeness to the gelation 
transition e is well below the Ginzburg-de Gennes crite- 
rion [24j which is a feature of the highest Mw polymers 
that distinguishes them from the others in the set. Hence 
the difference is likely to be due to non-mean-field behav- 
ior of these samples. Also, the four highest molar mass 
samples were prepared under slightly different conditions 
- where a partial reaction was carried out with stirring 
and, for the later part, the samples were reacted with- 
out stirring at a slightly elevated temperature Thus, 
our assumption of continuous stirred reaction may not 
be completely true for these samples. The dynamic ex- 
ponents calculated from our calculations match with the 
experimental findings in the relevant Mw range. 

To investigate the behavior of the dynamic exponents 
with average segmental lengths between branch points, 
we calculated the relaxation properties of a series of 
molecules in the ideal mean-field gelation ensemble. The 
dependence of the relaxation exponent u on N/N e falls 
about midway between the prediction of |17| and the phe- 
nomenological form of Q . We find that the viscosity ex- 
ponent becomes smaller as e is lowered. This is due to the 
dominance of supertube Rouse relaxation at long time 
scales for this class of polymers and, for small enough 
e, the viscosity exponent for any N/N e approaches the 
Rouse value applicable to the unentangled polymers. 

The recoverable compliance exponent t in our calcula- 
tions have values similar to those found in experiments. 
It is worth noting, however, that the magnitude of J°, 
when calculated from our algorithm, is found to be much 
larger than experimental values on similarly branched 
systems. Being the first moment of the relaxation mod- 
ulus G(t), the dominant contribution to J® comes from 
the long time behavior of G(t), so is very sensitive to 
the assumptions on which the relaxation dynamics of the 
very largest clusters in the ensemble is based. The com- 
putational scheme we used to follow the relaxation in 
the melt extrapolates ideas of dynamic dilation and su- 
pertube relaxation which originally were formulated for 
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linear or lightly branched systems to a highly branched 
system. In particular it assumes that the final supertube 
relaxation follows a Rouse scaling corresponding to a lin- 
ear object. The final relaxation, provided that the tail 
of the distribution is long enough, of largely unentan- 
glcd high molar mass molecules may find a faster route 
by showing a Zimm like relaxation, by which the largest 
clusters relax hydrodynamically in an effective solvent 
provided by the smaller clusters. In linear systems the 
transition molecular weight for this is the same as that 
for incomplete static screening of the larger molecules' 
self-interactions. Experimental results on model systems 
with high seniority and well characterized branching and 
molar mass are needed to quantitatively test the validity 
of the theory for accounting the long-time decay of stress 



in such highly branched systems. 

In summary, a numerical calculation of the entangled 
rheology of a series of mean-field gelation ensemble poly- 
mers provide a remarkable support of the accuracy of 
the hierarchical relaxation process suggested by the tube 
model. 
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